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Abstract 

We investigated the effect of deuteration on the vibrational ground state of the hydrated hydroxide 
anion using a nine-dimensional quantum dynamical model for the case of J = 0. The propagation 
of the nuclear wave function has been performed with the multi-configuration time-dependent 
Hartree method which yielded zero-point energies for the normal and fully deuterated species in 
quantitative agreement with previous diffusion Monte Carlo calculations. According to the zero- 
point energy the isotopomers having the hydrogen atom in the bridging position are more stable 
by about 1 kJ/mol as compared to the deuterium case. This holds irrespective of the deuteration 
state of the two OH groups. We also report the secondary geometric H/D isotope effect on the 
0-0 distance which amounts to an elongation of about 0.005 A for the symmetric isotopomers 
and 0.009 A in the asymmetric case. Finally, we explore the isotopomer sensitivity of the ground 
state tunneling splitting due to the torsional motion of the two OH groups. 
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I. INTRODUCTION 



Charged clusters with strong hydrogen-bonds have attracted considerable interest re- 
cently due to the progress of both, infrared spectroscopy and multidimensional quantum 
dynamical modeling l|, |2j. In particular the hydrated proton and its negative analogue, the 
hydrated hydroxide anion, have been investigated in quite some detail. For example, the 
Zundel cation, H 5 0^, became the paradigm for the predictive power of ab initio quantum 
dynamical calculations 3j, |J, |5J. The monohydrated hydroxide anion, H30^~, has been stud- 
ied by Johnson and coworkers using argon predissociation spectroscopy. In the range above 



3000 cm a sharp doublet has been detected 



at about 3650 cm which was subse- 



quently assigned to emerge from transitions between the HOOH torsion doublet and the 
"free" OH stretchings 8|. The spectral range from 1000-1900 cm -1 is dominated by a peak 
at 1090 cm -1 which was attributed to result from a combination of bridging hydrogen (BH) 
stretch, wag, and rock motions of the whole complex 9j. Effectively, this corresponds to a 
displacement of the BH away from the 0-0 axis. The region below 1000 cm -1 was finally 



addressed in Ref. 



101 ] . Besides two smaller features at 940 cm 1 and 995 cm 1 , with one 



being the second perpendicular mode of the BH, the spectrum is dominated by an intense 
and rather narrow peak at 697 cm -1 which has been assigned to the BH motion along the 
0-0 axis. 

The principal difficulty in assigning the H3O2 spectrum arises from the fact that similar 
to the Zundel cation one has to deal with a strong hydrogen bond in a floppy structure. 
In other words, the vibrational dynamics is rather anharmonic and any calculation bound 
to the harmonic approximation is likely to fail at least for the BH modes. Moreover, this 
system features a double minimum potential with a barrier of only about 75 cm -1 jll |. 
Therefore, the classical equilibrium structure will be nonsymmetric, e.g. [H-O-H- • ■ O- 
H]~, but as the zero-point energy (ZPE) including the BH motion is above the barrier 
the symmetric C2 structure H-O- • -H- • • O-H] - will be stabilized [ll, [l2, [3] as in other 
hydrogen-bonded clusters [iJ] . Of course, this affects especially the transition frequency of 
the BH stretching vibration drastically. Klopper et al. [ill] calculated a one-dimensional 
potential energy curve for this vibration which gave a ZPE about 350 cm -1 above the barrier 
and a fundamental transition frequency of about 1030 cm -1 which is lower than the harmonic 
prediction for a nonsymmetric structure (BH stretch coupled to the water bend around 



2 



1600-1700 cm -1 [l2j), but still far above the experimental value of 697 cm -1 lOj. This 
points to the need for a multidimensional treatment of the infrared spectrum as reported by 
J. Bowman and coworkers who used multimode reaction path vibrationally self-consistent 
: ield /configuration interaction (VSCF/CI) and diffusion Monte Carlo (DMC) techniques 
s[ Q]- F° r the specific case of the BH stretch these methods gave 741 cm -1 and 644 cm -1 , 
respectively. The accuracy (most modes within 38 cm _1 )of the VSCF/CI was recently 



confirmed by Lanczos diagonalization calculations in Ref. 

In this contribution we will focus on the effect of H/D isotopic substitution on the prop- 
erties of the hydrogen-bond in this complex. Our primary goal is to answer the question 
whether deuteration of the hydroxide ([D-O--- H- • • 0-H]~) or the "solvent" water ([H- 
O- ■ ■ D- ■ ■ 0-H]~ yields an energetically more stable structure. This issue is related to the 
problem of isotopic exchange equilibriums and associated fractionation factors which have 
been discussed for gas phase reactions of hydrogen-bonded ions [l?], [l8, Q as well as in 
aqueous 20|, l21| and other 22|, [23J solutions. An extensive account on the stability of various 
charge and neutral water clusters has been given by Scheiner and coworker employing the 
MP2 method together with a (small) 6-31+G** basis set 241 ] . The picture which emerged 
upon comparing harmonic ZPEs has been as follows: In neutral clusters like the water dimer 
the O- • • D- • • O hydrogen bond is stronger than the O- • • H- • • O one mostly due to the in- 
termolecular vibration of the bridging proton/deuteron. In charged clusters O- ■ ■ H- ■ ■ O 
wins as a consequence of the change in ZPE of the intramolecular water OH/D stretching 
vibration. In passing we note that for the Zundel cation this conclusion has been recently 
confirmed using anharmonic calculations^. 26]. For H^O^ the difference in ZPE is as small 
as about 50 cm -1 , although it was argued that this effect is enhanced by entropic effects. 
A related question could be asked for the doubly substituted case. Here, harmonic analysis 
predicts [D-O- • • H- • ■ O-D] - to be slightly favored over [H-O- • • D- • • O-D] - . It is important 
to emphasize that the conclusions for HaO^" have been drawn from harmonic calculations 
with respect to nonsymmetric equilibrium structure. In view of the discussion above this 
calls for having a second look at this problem from the perspective of a quantum mechanical 
treatment in full dimensionality. 

The change in strength of the hydrogen bond upon isotopic substitution is also reflected in 
geometric isotope effects (GIE) which is yet another manifestation of the multidimensional 



anharmonic nature of the potential energy surface (PES) 27, 



For weak hydrogen 
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bonds, having a double minimum PES, H/D substitution leads to a shorter 0-D distance as 
compared to 0-H. This in turn weakens the hydrogen bond thus causing the 0-0 distance 
to increase. For strong symmetric hydrogen bonds where the ZPE is above the barrier and 
the vibrational distribution has its maximum at the barrier top, H/D substitution reduces 
the width of this distribution which pulls the oxygens towards the deuterium, that is, the 
0-0 distance decreases. In the present paper we will address the GIE for all isotopomers 
and different degrees of deuteration on the basis of the full-dimensional ground state wave 
function, that is, at zero temperature. 

In the following Section III Al we will present a nine- dimensional (9D) Hamiltonian which 
describes the vibrational motion of H3O2 in full-dimensionality for total angular momentum 
equal to zero. This Hamiltonian is based on the PES developed by Bowman and coworkers 



15| • The related operator for the kinetic energy is given in the Appendix. The vibrational 



ground state is obtained by imaginary time-propagation using the Multi-configuration Time- 



dependent Hartree (MCTDH) method 29_|, |30[ ; numerical details of the calculation are given 



in Section III Bl In Section IIHI we will present results on the vibrational ground state of the 
different isotopomers and on the GIE. The paper is summarized in Section HVl 



II. THEORETICAL MODEL 



A. 9D Hamiltonian 



The choice of coordinates is crucial as it determines the strength of correlations between 
different degrees of freedom in the potential and kinetic energy operator. For the present 
case the issue is complicated as there is the possibility of large amplitude torsional motion. 



This made it necessary to use a reaction path approach in Ref. [151 ] . which combined the 
torsional reaction coordinate with orthogonal normal mode displacements taken with respect 
to the C2 transition structure. Here we will use internal coordinates, which is less restrictive 
but comes at the expense of a more complicated kinetic energy operator. After separating 
the total center of mass, the four Jacobi vectors shown in Fig. [1] will be used: Ri and 
R2 each connecting one oxygen and the "free" hydrogen atom, R4 connecting two centers 
of mass of the OH groups, and R3 connecting the shared hydrogen atom and the center 
of mass of the O2H2 fragment (accordingly for the deuterated cases). Assuming the total 
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angular momentum J 2 = the kinetic energy operator depends on the internal coordinates 
only. Based on these Jacobi coordinates the following nine internal coordinates are chosen: 
the lengths Ri, R 2 , and _R 4 , of the vectors Ri, R 2 , and -R 4 , respectively, the three Cartesian 
components of the BH, x, y, and z, being the components of R3, the angle 81 (8 2 ) between 
Ri and R$ (R 2 and R±), and the dihedral angle between the planes spanned by the 
vectors (Ri, R4) and (-R29 R4) which describes the torsional motion. Notice that due to 
its definition with respect to the Jacobi vectors, <p is slightly different from the torsional 



coordinate used in Ref . [15j . For the origin we use the center of mass and z is defined along 
the direction of i? 4 , i.e. it roughly corresponds to the BH stretch mode. 

In principle the exact kinetic energy operator can be obtained in terms of arbitrary 
coordinates by proper quantization of the classical Lagrangian. As it contains hundreds of 
terms one needs to compromise and we neglected the coupling between the angular momenta 
of the BH and the O2H2 fragment. For J 2 = this yields the kinetic energy operator given 
in the Appendix which contains 21 terms only. Notice that in the following numerical 
simulations we will use the new variables u% = cos 8i(i = 1, 2). 

The full- dimensional potential energy surface is constructed by cumulative expansion of 
different correlation orders [3l| . Following the strategy of Ref. |4J we will combine certain 
groups of coordinates to treat their correlation exactly. Specifically we have chosen the three 
groups g 1 = [Ri, R 2 , i? 4 ], g 2 = [ui, u 2 , (ft], and g 3 = [x, y, z]. The final expansion we used to 
generate the 9D PES reads 

V( gi ,g 2 ,g 3 ) = V + J2 y {1 \9i) + E ^tooty), « 

i i<j 

where gives the n-mode correlation between sets of coordinates. Notice, however, that 
this PES contains up to 6-mode correlations between individual coordinates. In Eq. (JT]) 
Vq is the energy of the reference geometry, g(°\ An obvious choice for this reference is the 
classical transition state geometry. We have performed a ground state calculation based on 
this reference and obtained an error of about 2-3 cm -1 per degree of freedom as compared 
to the Quantum Monte Carlo results for reported in Ref. [l5]. In a next step we 

have attempted to improve this result and found that one can get essentially a quantitative 
agreement if the quantum mechanical ground state expectation values (as obtained from the 
classical reference calculation) for the H 3 02~ coordinates are used to define the reference. 
For simplicity we have employed this reference for all isotopomers (see Table [T|) . 



For generating the PES in the given coordinates we have used the fitted CCSD(T)/aug- 
cc-pVTZ potential of Bowman and coworkers [15]. Two representative cuts of the PES are 
shown in Fig. [2] 



B. MCTDH Propagation 



The vibrational ground state of the different isotopomers has been obtained 
merical solution of the time- dependent Schrodinger equation in imaginary time 



the propagation we used the MCTDH method [29|, |30|] as implemented in the Heidelberg 



rom nu- 



32] . For 



program package 



331 ] . Here, the wavefunction is expanded as a superposition of Hartree- 



products consisting of time-dependent single particle functions (SPFs). A prerequisite for 
efficient propagation is a Hamiltonian which has product form. While the kinetic energy 
operator fulfills this requirement, the potential has been fitted to a sum of products using the 
POTFIT approach [30|. For the representation of the SPFs a discrete variable representation 
(DVR) has been utilized. For the torsional coordinate (f) an exponential DVR representation 
with periodic basis functions (eigenfunctions of Jr) was used. All other coordinates have 
been expressed via a harmonic oscillator DVR. The SPF basis functions cover the range 
energies below 10000 cm -1 ; the smallest natural orbital population in the ground state was 
0.0003. All parameters are compiled in Tab. [H 



III. RESULTS 



A. Vibrational Ground State 



Representative cuts of the full 9D ground state vibrational density are shown for the 
different isotopomers in Fig. [3] and the coordinate expectation values and their variances are 
compiled in Tabs. [Ill and IHII Fig. [3^ gives the densities along the shared proton coordinate 
z. For the symmetric cases [H-O- • • H- • • 0-H]~ and [H-O- • • D- ■ • O-H] - the distributions 
have their maxima at z = indicating the symmetrization due to ZPE which is above 
the classical barrier in Fig. [2h- As expected the width of the distribution narrows upon 
isotopic substitution. Inspecting Tab. [II] one finds, that this amounts to ~6% for the z 
coordinate, but to as much as 13 % for the bridging hydrogen's bending coordinates x and 
y. The distributions along the other coordinates are much less affected, e.g., the difference 
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for the torsion coordinate is not noticeable on the scale of Fig. [3b. Asymmetric isotopic 
substitution in [D-O- • • H- • ■ 0-H]~ introduces changes mostly at the substitution site, that 
is, R\ decreases and the associated width drops by 15%. Similar changes in the widths are 
observed for the angles U\ and (f) while the widths for the BH coordinates are essentially 
unaffected. The introduced asymmetry can also be deduced from the distribution along the 
coordinate in Fig. [2b. The widths for the bridging hydrogen do not change much in 
[D-O- • -H- • ■ 0-H]~ , however, as compared to [H-O- • • H- • ■ 0-H]~ the structure becomes 
asymmetric also insofar as the hydrogen atom moves toward the OH group and out of the 
0-0 axis, i.e. [D-O- • -H-O-H] - is formed. Related trends are observed if one starts from 
[D-O- • • D- ■ • 0-D]~ and makes D/H substitution. The definition of the coordinates depends 
on the isotopomer which affects the asymmetric cases. Therefore, we also give the values 
of internal coordinates in the table captions in order to support the conclusion mentioned 
above. 

Next we focus on the ZPEs of the different isotopomers which are analyzed in terms of 
the contributions of the different parts of the Hamiltonian in Tab. IIV1 First we notice, that 
the ZPE of [H-O- • • H- • • 0-H]~ is calculated as 6606 cm -1 which essentially reproduces the 
DMC result (6605± 5 cm -1 ) obtained on the fully coupled PES [l5j. The set g 1 gives the 
largest contribution to the ZPE as it contains the high frequency OH stretching vibrations. 
Further, the strongest correlations between different sets of coordinates are those involving 
the BH motion. For the fully deuterated case we obtain 4481 cm -1 which is also in accord 
with the DMC result (448 7± 5 cm" 1 ). 

Next we discuss the general trend in ZPE change upon H/D substitution. Inspecting 
Table IIVI we observe that replacing H by D in one of the OH groups lowers the ZPE by 
about 600 cm -1 irrespective whether the other OH group is deuterated or not. Replacing H 
by D in the bridging site lowers the total ZPE by about 520 cm -1 only, again irrespective of 
the OH groups. In other words, in terms of the ZPE the H-bond is about 80 cm -1 stronger 
than the D-bond. Let us have a more detailed look at the single substitution case. Here, we 
find that in terms of ZPE [D-O- • • H- • • O-H] - is more stable than [H-O- • • D- • • O-H] - by 
82 cm -1 . First of all, this confirms the qualitative result of Scheiner and coworker jiyj] who 
obtained a value of 52 cm -1 based on the harmonic approximation of the PES around a 
nonsymmetric structure. Their normal mode treatment led to the conclusion that it is the 
intramolecular water OH(D) stretching which is responsible for the increased stability of [D- 
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O- • ■ H- ■ • 0-H]~ . However, as discussed before the structure is symmetrized by ZPE and 
its analysis in terms of a water molecule hydrogen-bonded to a OH(D) _ is not adequate. 
Instead inspecting Tab. IIVI we find the following behavior: First, ZPE changes in [D- 
O- • ■ H- ■ ■ 0-H]~ are due to the set g l , i.e. mostly the OH-stretching at the deuteration 
site, while in [H-O- • • D- • • 0-H]~ it is set g 3 involving the BH. We also notice that the net 
effect of all single mode potentials is 37 cm -1 by which the D-bond would be more stable 
whereas the correlation energy difference of -119 cm' 1 finally leads to the preference for 
the H-bond. Although we should emphasize that the actual values of the different energy 
contributions depend, of course, on the chosen reference geometry for the PES expansion, it 
seems plausible that, for instance, the correlation between the bridging hydrogen set g 3 and 
both sets g l and g 2 plays a significant role for the ground state wave function. As a caveat 
we note that these findings do not necessarily imply a simultaneous coupling between all 
modes of the given sets. 



B. Secondary geometric isotope effect 

Having at hand the ground state wave functions we can calculate the secondary GIE, i.e., 
the change of 0-0 distance upon isotopic substitution. Since does not fully correspond to 
the 0-0 distance, Ro-o, (cf. Fig. [Q the respective operator had to be expressed in terms 
of our model coordinates. The resulting expectation values for the different isotopomers 
are compiled in Tab. |Vj The classical value for Ro-o m HOHOH at the MP2/aug-cc- 
pVTZ transition state is 2.446 A. In the quantum case this value increases to 2.492 A as 
a consequence of zero-point vibration. This value is reduced to 2.488 A for substitution of 
the BH, i.e. in the HODOH case. In all cases the substitution of the central BH leads to a 
reduction of the 0-0 distance. For the symmetric cases the reduction is about 0.005 A while 
for the asymmetric case it is 0.009 A. Such a bond compression due to reduced zero-point 
vibration, i.e., localization of the wave function, is typical for strong hydrogen bonds. 

The changes reported here are in good agreement with the DMC calculations in Ref. 
for H 3 2 and D^O^ , only the absolute values of the DMC Ro-o are larger by 0.005 
A. In passing we note that the fully deuterated case has also been investigated using the 



finite temperature path integral method at a lower level of quantum chemistry 34| . For this 
situation Tachikawa and coworker obtained 2.498 A for H3O2 and 2.504 A for D 3 0^~, that 
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is, the opposite trend which has been explained by the bimodal character of the calculated 
distribution for at the given level of theory. 

IV. SUMMARY 

The monohydrated hydroxide anion is a prototype for a strongly hydrogen-bonded low- 
barrier system whose structure is symmetrized by zero-point vibrational motion. An ac- 
curate theoretical prediction of the associated ground state wave function requires to treat 
the dynamical problem in full-dimensionality. Using the CCSD(T) potential energy surface 
of Bowman and coworkers [15] we have shown that the MCTDH approach for wave packet 
propagation [29J can meet this challenge. Comparing the vibrational ground state wavefunc- 
tions and ZPEs for different isotopomers the following conclusions could be drawn: First, 
in accord with the general view that for ions bridging donor and acceptor by a hydrogen 
atom is more favorable than bridging by a deuterium. The general trend for the considered 
systems is that in terms of ZPE the H-bonds are about 80 cm^ 1 (~1 kJ/mol) more favorable 
than the D-bonds, irrespective the deuteration state of the OH groups. Specifically, we find 
[D-O- • -H- - ■ 0-H]~ to be energetically more stable than [H-O- • • D- - • 0-H]~ by 82 cm -1 . 
Although this seems merely to confirm the results of the harmonic analysis reported in Ref. 



24j , the present full- dimensional treatment is providing a more realistic physical picture by 
accounting not only for the symmetrization of the structure due to zero point motion but 
also for the anharmonicity of the potential energy surface. Viewed from the perspective of 
harmonic vibrations it is not only the loss of an intramolecular OH vibration of the water 
molecule which reduces the ZPE more than the loss of the intermolecular bridging hydrogen 
vibration. In fact we find that correlations between the bridging H/D atom motion and 
especially the 0-0 stretching contribute significantly to the ZPE. 

The second result concerns the H/D isotope effect on the heavy atom 0-0 distance. 
Here, we find that as compared with the classical prediction, zero-point vibrational motion 
increases -Ro-o m the H-bonded isotopomers. Due to the reduced wave function derealiza- 
tion in the D-bonded forms Rq-q decreases again, although not below the classical value. 
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Appendix 

The kinetic energy operator is derived by direct quantization of the classical kinetic energy 
Tdass = MQ, where Q is a vector including all generalized coordinates and f means 
Hermitian conjugate which is equivalent to transpose in the classical kinetic energy. The 
corresponding momentum operator is defined by P = —ihd/dQ and hence the quantum 
kinetic energy operator is T = |pt M X P . The general method to derive the Hermitian 
conjugate of the momentum will be detailed elsewhere. The classical kinetic energy in 
laboratory reference frame is obtained by conventional mechanics as a sum of two parts 
Tciass = T<Aaa&{Rz) + Tciass(-Ri) R2, R4) ■ For the latter part in the sum we have ignored the 
global rotation of the O2H2 fragments. In other words the angular momentum coupling 
between the central hydrogen and the O2H2 fragments has been neglected in the expression 
for the exact kinetic energy operator for total angular momentum J 2 = 0. For the simulation 
an additional normalization transform was performed to slightly reduce the numerical effort. 

This gives the following 9D kinetic energy operator: 

T = T 1 + T 2 + T 3 

with 

h 2 d 2 h 2 d 2 

1 ~ 2fi 1 dR 2 2fi 2 dR 2 2fi 4 dR 2 
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, 1 1 \ d „ 2 , d ( 1 1 \ d 2N d 
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2^1 ^T^u 2 \d(j) "V 'duz du 2 

h 2 d 2 h 2 d 2 h 2 d 2 



2/i 3 dx 2 2/jl 3 dy 2 2// 3 dz 2 

The non-Euclidean normalization according to the volume element is d T = 
dRidR2dR 4l dxdydzduidu2d(f). The reduced masses are defined as follows: HOHOH 
Hx = /i 2 = m R m / (m H + m ), A*3 = 2m H (m H + m )/ (3m H + 2m ), /i 4 = (™ H + m )/2; 
HODOH - /i 3 = 2mD(mn + m )/(m D + 2m H + 2m ); HOHOD - - Hi = m D m / (m H + mo) 
and /X3 and /x 4 change correspondingly. For DODOD - the corresponding masses of HOHOH 
are modified by replacing mn by m^. In the same way we get the masses for DOHOD - 
by exchanging mn and in HODOH~, and similarly one can obtain HODOD from 
DOHOH . 
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TABLE I: MCTDH parameters for the imaginary time propagation (lengths in A, A^dvr : number 
of DVR points; Wspf : number of SPFs.) and reference geometry, g'°J , for the PES expansion in Eq. 
([1]) (X=H/D). Note that the coordinates of the reference geometry are based on expectation values 
for the HOHOH case. The actual different values reported here for the different isotopomers are 
due to the fact that the Jacobi vectors defining the coordinate system are mass-dependent. The 
classical transition state (TS) reference for HOXOH - is given as well. 





i?l i?2 R& 




x y z 


min. grid 
max. grid 
-Wdvr 


0.740 0.740 2.275 
1.322 1.322 3.069 
11 11 16 


-1.0 -0.7 
0.7 1.0 2vr 
13 13 17 


-0.635 -0.635 -0.529 
0.635 0.635 0.529 
13 13 17 


iVsPF 


5 


10 


13 


g®: HOXOH" 
flf<°): DOXOD" 
gf(°): DOXOH" 
TS: HOXOH 


0.979 0.979 2.519 
0.979 0.979 2.547 
0.979 0.979 2.532 
0.962 0.962 2.477 


-0.25 0.25 7T 
-0.288 0.288 vr 
-0.269 0.269 vr 
-0.282 0.282 vr 


0.0 0.0 0.0 
0.0 0.0 0.0 
-0.012 0.021 0.043 
0.0 0.0 0.0 
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TABLE II: Coordinate expectation values and their variances (in A) for the 9D ground states 
of H3O2 and different single D substitution isotopomers. The expectation values can be trans- 
lated into internal coordinates which gives for DOH*0*H- : i? O -H*=1.270 A, # *-H*=l-226 A, 
angleo-H* -o* = 178.8° . 





HOHOH 


HODOH 


DOHOH 




mean 


variance 


mean 


variance 


mean 


variance 


JrL\ 


0.981 


0.070 


n no 1 

0.981 


0.070 


0.97b 


0.059 




0.981 


0.070 


0.981 


0.070 


0.981 


0.070 


R 4 


2.521 


0.063 


2.515 


0.064 


2.536 


0.065 


Ul 


-0.259 


0.173 


-0.261 


0.172 


-0.267 


0.149 


U2 


0.259 


0.173 


0.261 


0.172 


0.269 


0.174 


<t> 


3.138 


1.240 


3.136 


1.238 


2.886 


1.258 


X 


0.000 


0.125 


0.000 


0.107 


-0.022 


0.126 


y 


0.000 


0.120 


0.000 


0.104 


0.029 


0.123 


z 


0.000 


0.150 


0.000 


0.140 


0.029 


0.152 
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TABLE III: Coordinate expectation values and their variances (in A) for the 9D ground states 
of and different single H substitution isotopomers. The expectation values can be trans- 

lated into internal coordinates which gives for DOD*0*H~: i?o-D*=l-263 A, i?o*-D*=l-225 A, 
angleo-D*-0* =178.9°. 





DODOD 


DOHOD 


DODOH 




mean 


variance 


mean 


variance 


mean 


variance 


JrL\ 


0.976 


U.U59 


0.976 


U.U59 


0.97b 


0.059 




0.976 


0.059 


0.976 


0.059 


0.981 


0.070 


R 4 


2.544 


0.064 


2.549 


0.064 


2.528 


0.064 


Ul 


-0.281 


0.148 


-0.279 


0.148 


-0.270 


0.146 


U2 


-0.281 


0.148 


0.279 


0.148 


0.272 


0.172 


<t> 


3.139 


1.290 


3.135 


1.293 


2.889 


1.254 


X 


0.000 


0.120 


0.000 


0.135 


-0.021 


0.109 


y 


0.000 


0.104 


0.000 


0.120 


0.028 


0.105 


z 


0.000 


0.142 


0.000 


0.151 


-0.027 


0.140 
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TABLE IV: Energy expectation values (in cm' 1 ) of the different isotopomers in the vibrational 
ground state. The subscripts refer to single sets for one set operators or to pairs of sets for the 
two-set operators, e.g.. Vu = (V^ (g 1 , g 2 ))- Note that Vq=462 cm -1 in Eq. (pQ). 





T To To T 
J-1 J 2 ^3 1 


Vx V 2 


V 3 


vi2 




Vnn V 

V23 V 


T-i-V 


HOHOH" 


1995 338 935 3268 


1976 425 


944 


-5 


-80 


-384 3338 


6606 


HODOH" 


1997 349 646 2992 


1964 423 


681 


-5 


-71 


-359 3095 


6087 


DOHOH" 


1741 299 940 2980 


1731 392 


994 


-4 


-122 


-428 3025 


6005 


DODOH" 


1745 309 651 2705 


1698 374 


685 


-5 


-68 


-368 2778 


5483 


DOHOD 


1489 259 941 2689 


1459 385 1005 


-6 


-93 


-496 2716 


5405 


DODOD 


1491 268 653 2412 


1445 378 


740 


-6 


-83 


-467 2469 


4881 
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TABLE V: Expectation values for the 0-0 distance (in A) and reductions for substitution of BH 
for the different isotopomers. The classical value at the transition state is 2.446A. 



case (Ro-o) 


case (Ro~o) 


case (Ro-o) 


HOHOH 2.492 
HODOH 2.488 


DOHOH 2.495 
DODOH" 2.486 


DOHOD 2.493 
DODOD 2.488 
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FIG. 1: Classical equilibrium structure of the H3O2 anion and the four Jacobi vectors used for 
defining the nine internal coordinates (see text). For the PES expansion in Eq. (pQ) we have used a 
reference geometry based on the coordinate expectation values as obtained for H3O2 (cf. Tab. HJ). 
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FIG. 2: (a) PES along 0-0 distance (-R4) and proton transfer coordinate z (contour lines at (in 
cm" 1 ): 230, 260, 300, 600, 900, 1300, 1700, 2100, 2500, 3000). (b) Potential energy curve along 
the coordinate <j) which for the symmetric case corresponds to the usual torsion angle. In both 
panels all other coordinates have been kept frozen at the transition state geometry. 
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Torsion angle (degree) 



FIG. 3: Reduced probability density for the vibrational ground states of the different isotopomers 
for the proton transfer (a) and the torsional coordinate (b). 
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